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The two-dimensional Ising model on a distorted Kagome lattice is studied by means of exact solutions and the 
tensor renormalisation group (TRG) method. The zero-field phase diagrams are obtained, where three phases 
such as ferromagnetic, ferrimagnetic and paramagnetic phases, along with the second-order phase transitions, 
have been identified. The TRG results are quite accurate and reliable in comparison to the exact solutions. In a 
magnetic field, the magnetization (m), susceptibility and specific heat are studied by the TRG algorithm, where 
the m = 1/3 plateaux are observed in the magnetization curves for some couplings. The experimental data of 
susceptibility for the complex Co(N 3 ) 2 (bpg)- DMF 4/3 are fitted with the TRG results, giving the couplings of 
the complex J = 22K and J' = 33K. 
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I. INTRODUCTION 

Kagome lattice, one of the most interesting frustrated spin 
lattices, has attracted much attention both experimentally and 
theoretically in recent years. The Heisenberg model on a 
Kagome lattice is likely a candidate for finding a spin liquid 
state in two-dimensional (2D) spin systems. Some numerical 
works have revealed that the system possesses a magnetic dis- 
ordered ground state. 1 ' 2 Nevertheless, the nature of its ground 
state is still an open question. 3 Recently, a number of spin 
systems, such as volborthite Cu 3 V 2 07(OH) 2 -2H 2 (Refs. 
043), [H3N(CH 2 ) 2 NH 2 (CH 2 ) 2 (NH3] 4 [FeJ ) I F ls (S04)6]-9H 2 
(Ref. S) and Co(N 3 ) 2 (bpg> DMF 4/3 (Ref. @), are found to 
form a Kagome lattice with distortions, where the structural 
distortions give rise to two different exchange couplings J and 
J'. Such a spatially bond anisotropic spin lattice can be called 
a distorted Kagome (DK) lattice, as schematically depicted in 
Fig. Q] A distortion-induced magnetization step at small fields 
and the 1/3 magnetization plateau on the DK lattice have been 
observed in experiments. 7 In order to explain the experimen- 
tal results, the quantum and classical Heisenberg models on 
the DK lattice are considered42iii In most cases, the interac- 
tions between spins are usually of Heisenberg type on a DK 
lattice 10 . However, in some materials the spin-spin couplings 
are anisotropic, and even of Ising-type in particular situations. 
For instance, in complex Co(N 3 ) 2 (bpg)- DMF4/ 3 , Co +2 ions 
form a spin- 1/2 DK lattice and might be coupled by Ising-type 
interactions at low temperatures. 12 Therefore, it should also be 
necessary to pay more attention on the 2D Ising model with 
the DK lattice, especially when an external magnetic field is 
present, where the works in the literature are still sparse. 

In this article, we shall focus on the Ising model on the 
DK lattice, where the thermodynamics and magnetic proper- 
ties will be carefully studied by means of exact solutions and 
a numerical method. In the present model, the quantum fluc- 
tuations are completely suppressed, and only the thermal fluc- 
tuations are considered. It should be stressed that the present 
model for h — can be exactly solved, but for h + the 




FIG. 1: (Color online) (a) The distorted Kagome lattice, where J 
and J' denote the different nearest neighbor couplings. The dots in 
the center of small triangles and lines connecting them form a tensor 
network presented by Eq. l[3}. (b) and (c) show two degenerate spin 
configurations of the ferrimagnetic structure, where up and down ar- 
rows represent the spin-up and spin-down states, respectively. 



exact solution is not available at present and the numerical 
method should be involved in. For this reason, a recently de- 
veloped tensor renormalization group (TRG) method i 3 ^— is 
employed to investigate the thermodynamic properties of the 
system. The cases with different couplings J and J' and both 
with and without a magnetic field h will be discussed. The 
TRG results, being consistent with the exact solutions, reveal 
that the system has ferromagnetic, ferrimagnetic, and para- 
magnetic phases in the phase diagram, where the paramag- 
netic phase can exist at T — in the strongly frustrated re- 
gion. Moreover, the magnetic order-disorder phase transitions 
separating these three phases are disclosed. In ferrimagnetic 
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FIG. 2: (Color online) The zero-field phase diagrams for different 
couplings: (a) J > and (b) J < 0. The phase boundaries are deter- 
mined by Eq. (|2j. In the phase diagrams, FM means ferromagnetic 
phase, FI represents ferrimagnetic phase, and PM is paramagnetic 
phase. 



the free energy per triangle in the thermodynamic limit as 
T 



f 



t \ d0 < 



ln[a + 2b cos(0) 



16tt 2 J_„ 
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where 

2 2 2 ? 

a = Wj + W 2 + w 3 + w 4' 

= (jJ\U)t, - U>2U>4, 

C — Cl)\U>4 — LL>2ti)3, 

d — W3W4 - 

e — 0)30)4 — u>5tO(>, 

u>i = 2exp(-27/r)[l +exp(2y/r)cosh(27'/r)] 2 , 

a>2 = oji - 8 cosh(2J'/T), 

W3 = a>4 — W5 = = exp(27/r) cosh(47'/r) - exp(2//T), 

w 7 = w 8 = wi - 4exp(-2//T), 



and paramagnetic phases, the 1/3 magnetization plateaux are 
seen. A magnetization step at an infinitesimal field is observed 
for the paramagnetic phase at low temperature. The specific 
heat no longer possesses a divergent peak once the external 
field is switched on, implying the absence of phase transi- 
tions at h + 0. In addition, we shall also make an attempt 
to fit the experimental data of susceptibility for the complex 
Co(N3)2(bpg)- DMF4/3 with the TRG results so as to estimate 
the exchange couplings in this complex. 

The other parts of this article are organized as follows. In 
Sec. II, the exact solutions and phase diagrams are presented. 
The TRG method is introduced in Sec. III. The specific heat 
without a magnetic field is explored in Sec. IV. Magnetization, 
susceptibility, and a comparison to the experimental data are 
shown in Sec. V. Sec. VI contains the results of specific heat 
in an external magnetic field. The summary and discussions 
are given finally. 



II. EXACT SOLUTIONS IN ZERO MAGNETIC FIELD 

The Hamiltonian of the system under interest has the form 

of 

H = r 2 SiSj+s X SiS ' +J X SiSj-hJ^St, 

<ieajefi> <ieajey> <iefij£y> i 

(i) 

where 5, is an Ising spin with two discrete values ±1, and the 
whole lattice can be divided into three sublattices labeled as a, 
(3, and y. The spin-spin coupling terms are restricted to nearest 
neighbor sites, J and J' are different nearest neighboring cou- 
plings as shown in Fig. \T\ where J(J') > and < represent 
antiferromagnetic and ferromagnetic couplings, respectively, 
h is the uniform external magnetic field, and gfig — 1 is as- 
sumed. 

Let us first utilize the exact mapping of 2D Ising model onto 
the 16- vertex model to present the exact solutions on the DK 
lattice. Following Refs. Il6l - [l8l for h — 0, we can write down 



where k B — 1 is presumed. It is straightforward to readily ver- 
ify that these w's satisfy the free-fermion conditions, 16 show- 
ing that the Ising model on the DK lattice is exactly solvable. 
The critical temperature T c at which the phase transition takes 
place is determined by the following equation 

1 + 4 exp(2//T c ) cosh(27'/r c ) - cosh(47'/r f ) = 0. (2) 

Two different cases are presented in Figs. [2] (a) and (b). 
When J > 0, the magnetic ordered phases appear when 
\J'/J\ > 1, which can be recognized as a ferrimagnetic phase 
for J' I J > 1, and a ferromagnetic phase for J' I J < -1. This 
is obtained by checking the magnitude of spontaneous mag- 
netization [see Figj3] (c) below]. The disordered paramag- 
netic phase separates the two ordered phases, where the phase 
boundaries are determined by Eq. (0. For a small T c , Eq. (O 
can lead to a simple expression T c /J m 2(\J'/J\ - 1)/ ln4, that 
gives the straight phase boundaries in Fig. [2] (a). Notice that 
the paramagnetic phase exists even at zero temperature owing 
to the frustration. When J > and \J'/J\ e [0,1], the spin 
surrounded by J' couplings on a sublattice is free to flop up 
or down without costing energy. Meanwhile, the ground-state 
spin configurations on fi and y sublattices are highly degener- 
ate. Hence, the total degeneracy is K — 2 N ^ i+Nc , with the 
number of total sites, and N c the number of chains consisting 
of spins on f3 and y sublattices (the horizontal lines in Fig. Q]). 
This superdegenerate state at zero temperature connects con- 
tinuously with the paramagnetic phase at finite temperatures. 
No phase transition appears, and the system is disordered at 
all temperatures. This observation can also be manifested in 
Fig. [3] (a), where no singularities exist in the specific heat for 
\J'/J\ = 0.4,0.8, 1.0. However, for \J'/J\ > 1, as Fig. 0(a) 
indicates, there exists a ferromagnetic (J' < 0) or ferrimag- 
netic (J' > 0) state at low temperatures, and these ordered 
phases would be destroyed through an order-disorder phase 
transition with increasing temperature. Correspondingly, the 
specific heat for \ J'/J\ = 1.2 in Fig. [3] (a) shows a divergent 
peak, which is a typical character of second-order phase tran- 
sition. When J < [Fig. 12(b)], the situations are similar, 
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FIG. 3: (Color online) Temperature dependence of specific heat and 
magnetization for different coupling ratio J' I J at h = 0. The TRG re- 
sults (open symbols) along with the exact solutions (solid and dashed 
lines) are presented for (a) J > and (b) / < 0. In (c) and (d), the 
magnetization m is plotted for J > and J < 0, respectively, where 
a magnetic order-disorder transition is clearly seen. The insets illus- 
trate the critical behaviors of m near T c . 

and the system is ordered (ferromagnetic or ferrimagnetic) at 
low temperatures except for the case J = -1,7' = 0, where 
the model is degenerated into the decoupled one-dimensional 
Ising chains, which has T c — and thus is disordered at any 
finite temperature. 

HI. TRG ALGORITHM 

Exact solutions can offer us a reliable phase diagram of the 
model. However, some other quantities such as the magnetiza- 
tion m and specific heat in nonzero magnetic fields cannot be 
obtained within the above framework. Hence, we adopt the re- 
cently proposed TRG numerical algorithm. The TRG method 
is first introduced to calculate the 2D classical modelS j 13 i 19 and 
then generalized to study 2D quantum spin modelsi 14,15 ' 20 i 21 
The principal idea of TRG algorithm is to express the partition 
function (or the expectation value of quantum operators) as a 
tensor network, and then utilizes the coarse-graining and deci- 
mation procedures to approximately obtain the results. TRG is 
an efficient method both for classical and quantum spin mod- 
els. 

The first step is to replace each triangle on Kagome lat- 
tice by a tensor, as shown in Fig. Q] (a). The energy of 
each triangle in an external magnetic field h is e A (s\, S2, S3) = 
J's t S2 + J's\Sj, + JS2S3 - hh(si + S2 + S3). We introduce 
a three-order tensor T^}^ S} = wp{—s A {s\,S2,s%)IT), where 
A(B) means down (up)-pointing triangle in Fig. [T|(a). These 
tensors form a honeycomb lattice, and the partition function 
can be expressed as 

Z= 2j exp{-[e A (*i,5 2 ,*3) + £aOi,*4,*5) + -]/T] 

Sl,S2,Ss,... =— 1,1 

E T L^ T l = tTHT A T B . ..), (3) 

Si,S2,Sy,...=— 1,1 




0.0 0.5 1.0 1.5 



FIG. 4: (Color online) The specific heat and magnetization m as func- 
tions of temperature with \J' /J\ = 0.04 at h = 0. (a) Specific heat, 
where three peaks appear, one of which is divergent; (b) Magnetiza- 
tion m and sublattice magnetization m a , where \m a \ decreases rapidly 
(but does not vanish) around the peak position of the specific heat 
at low temperature, which is also revealed as a local maximum of 
dm a /dT in the inset. 



where tTr represents the tensor trace. Hence the problem of 
solving the partition function of Ising model on a DK lat- 
tice is equivalently transformed into a honeycomb tensor net- 
work problem, which can be efficiently evaluated through the 
rewiring and coarse-graining iterations (see the details in Ref. 
fl4l) . Upon obtaining the partition function Z, other thermo- 
dynamic quantities can be evaluated straightforwardly. Al- 
ternatively, we can also introduce some impurity tensors in 
the tensor networks to achieve this goal. For example, in 
order to calculate the magnetization m, an impurity tensor 
T 1 "^ = ( £i± f ±£i )exp[-l/r £A (s 1 , j 2 , s 3 )] can be introduced. 
By replacing one tensor T A ^ B in Eq. (0, we can get the mag- 
netization per site 

tTr(T ,m T B T A ...) 

m = -. (A) 

Z 

In the following, the second scheme is adopted for evaluating 
the thermodynamical quantities, such as the magnetization m, 
energy per site e, etc. 

In our calculations, the number of coarse-graining iterations 
is generally taken as 20, i.e., the total sites of DK lattice under 
investigation is 3 22 ~ 3 X 10 10 , which is close to the thermody- 
namic limit. In addition, the periodic boundary conditions are 
adopted during the simulations. The initial bond dimension D 
of tensor T is chosen as 2 owing to the two states (spin-up and 
-down) of the Ising spins. With the coarse-graining procedure, 
the bond dimension will increase dramatically, and hence we 
have to make a truncation and reserve a finite dimension D c . 
In our calculations, D c is taken as high as 18, and the conver- 
gence with various D c has always been checked. 
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IV. SPECIFIC HEAT AND PHASE TRANSITIONS 

When h — 0, both exact solution and TRG method can be 
utilized to evaluate the specific heat. In Figs. [3] (a) and (b), 
the TRG results are plotted by symbols, while the exact re- 
sults by lines. Excellent agreement can be observed, except 
for the region around the critical point where a divergent peak 
occurs. Another character is that the specific heat at zero field 
is independent of the sign of coupling but is relevant to 
the magnitude of |/'|. In Fig. 0(a), when \J'\ is small, there 
is only one round peak in the specific heat. By tuning |/'| to 
approach J from below (e.g. \J'\/J - 0.8), a new round peak 
appears at low temperature, which disappears when |/'| = J 
and the system again exhibits a single round peak. These ob- 
servations imply that there exist no phase transitions when 
|/'| < /, which is owing to the strong frustration, and is in 
accordance with the phase diagram in Fig. [2] (a). Further- 
more, if \J'\ exceeds J [as \J'/J\ = 1.2 in the Fig. [3] (a)], 
a divergent peak emerges, implying the occurrence of phase 
transition. In Fig. [3(b), a typical curve of specific heat with 
\J' /J\ = 1 is shown. A divergent peak occurs at the transi- 
tion temperature T c « 2.14, again in agreement with the ex- 
act solution (T c = 4/ln(3 + 2 V3)). Moreover, the specific 
heat is logarithmically divergent at the critical point because 
of 2 exp(2//T c ) + cosh(47'/r e ) - 1 + 0, as shown in Figs. [3] 
(a) and (b). 

The phase transitions can also be verified by studying the 
order parameter, i.e., the magnetization per site m defined in 
Eq. ©. In Figs. [3] (c) and (d), when J = 1, J' /J > 1 
or J = -1,7' > 0, m = 1/3 at T = and remains finite 
at small temperatures. This nonzero spontaneous magneti- 
zation implies the existence of a ferrimagnetic phase; while 
J = 1, J' I J < — 1 or / = -I, J' < 0, the magnetization 
starts from m = 1, and the system is in a ferromagnetic phase 
when T is smaller than the critical temperature T c . With in- 
creasing temperature, m decreases steeply to zero in the vicin- 
ity of T c , showing a order-disorder phase transition happens. 
In addition, the critical behavior of m near T c has been in- 
vestigated. In the insets of Figs. [3] (c) and (d), in terms of 
m cc (T c - T) n , the fittings in different cases coincidentally 
give r] - 1/8, which is the same as that of Ising model on a 
square lattice. 22 The phase transition occurring at T c probably 
falls into the universality class of 2D Ising models. 

Another interesting case is shown in Fig. |4](a), where the 
temperature dependence of the specific heat is presented for 
J < and <K \J\. One may note that there are three peaks, 
including two round peaks and a divergent one. Both exact 
solutions and TRG method give the same results. In order to 
investigate the origin of each peak, the TRG method is uti- 
lized to calculate the magnetization m and sublattice magne- 
tization m a . As shown in Fig. [4] (b), m behaves rather dif- 
ferently for /' = 0.04 and -0.04, although the specific heat 
coincides for both. When /' = -0.04, the ground state is fer- 
romagnetic, and m decreases monotonously with increasing 
temperature and vanishes sharply at critical temperature T c . 
The case with J' = 0.04 is more interesting, where the system 
possesses a ferrimagnetic ground state with m — 1 /3 at T — 0. 
With increasing temperature, m first increases until the tem- 
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FIG. 5: (Color online) In (a) and (b), the magnetic curves for differ- 
ent couplings for J > and J < are shown, respectively, where 
\T/J\ = 0.2 and D r =18. Inset in (b) presents the magnetic curves 
with different temperatures below and above T c . In (c) and (d), the 
ground-state phase diagrams on J' - h plane are presented. 



perature is close to the critical point T c , and then goes down 
steeply to zero. In order to understand this peculiar behavior, 
we have plotted the sublattice magnetization m a as a func- 
tion of T in Fig. H](b). In the ferrimagnetic case, m a is aligned 
anti-parallel with the spins on the other two sublattices, and its 
magnitude decreases rapidly with increasing temperature ow- 
ing to the coupling J' weak. Hence, m = (-\m a \ + mg + m r )/3 
would first increase until the temperature approaches to T c , 
where m a , nip, and m y disappear simultaneously. Moreover, 
as the inset shows, the first-order derivative dm a /dT has a 
round peak at the temperature T r , which coincides with the 
low temperature peak position of specific heat. Although m a 
decreases rapidly around T r , it does not vanish. It should be 
pointed out that T r is not a critical point, as the specific heat 
shows only a round peak and never diverges at T r . 



V. MAGNETIZATION AND SUSCEPTIBILITY 
A. Magnetization Plateaux and Ground State Phase Diagrams 

When the external magnetic field is switched on, the exact 
solution in Sec. II no longer works. The TRG method, which 
has been verified to be accurate and reliable in the previous 
sections, is utilized to study the response of the system to an 
external magnetic field. 

Let us first focus on the magnetization, where the 1 /3 mag- 
netization plateaux are obtained, as shown in Fig. [5] When 
J > and J' I J e [-1,1] in Fig. |5] (a), the system is in 
a paramagnetic phase at all temperatures. An infinitesimal 
small magnetic field can polarize the free spin on a sublattice 
at T = 0, and hence drives the ground state to a ferrimag- 
netic state with m = 1/3 after a magnetization jump, and the 
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FIG. 6: (Color online) Temperature dependence of zero-field suscep- 
tibility for different J and /'. (a) x diverges at T = and^-r con- 
verges to 1/3 as the inset shows; (b) x shows divergent peaks at the 
critical temperature T c , where the magnitude of x with J = J' = -1 
has been divided by two. The susceptibility is calculated at h = 0.01 
and, the convergence with various small fields has been checked. 



1/3 plateaux appear in the magnetization curves. At finite but 
small temperatures, these plateaux are still present. This field- 
induced 1/3 plateau ferrimagnetic phase is highly degenerate, 
and the degeneracy is K — 2 Nc , where N c is the number of 
independent spin chains in the system. One of the degener- 
ate spin configurations is shown in FigJT] (b). When the field 
is larger than a critical field h c , the spins on /3 and y sublat- 
tices align parallel instead of antiparallel, and the system has 
the saturated magnetization m = 1, leading to a ferromagnetic 
spin configuration. The energy difference per site between the 
polarized ferromagnetic and the plateau ferrimagnetic state is 
5e = i(J' + J). The Zeeman energy 6e z = at the criti- 
cal magnetic field h c has to compensate this energy difference. 
5e z + be — leads to h c = 2(7 + 7'), as verified in Fig. |5](a), 
where the critical field h c increases with enhancing the cou- 
pling J'. In addition, there exists a notable difference between 
the magnetization curves of 7' — 1.2 and others in Fig. |5](a). 
The former starts from a nonzero spontaneous magnetization 
(m = 1/3) owing to its ferrimagnetic ground state instead of 
a paramagnetic one. The ground state spin configuration for 
7' > 1 is illustrated in Fig. [T|(c), which is ferrimagnetically 
ordered. Considering the spontaneously broken Z2 symmetry, 
this ground state is no longer degenerate. Hence, the width of 
1/3 plateau does not obey the relation mentioned above, but 
has another relation in the ferrimagnetic case to be discussed 
below. 

For 7 < 0, there exist ferromagnetic (7' < 0) and ferrimag- 
netic (/' > 0) ground states. The cases with 7' > possess 
m — 1/3 plateaux, as seen in Fig. 0(b). Comparing with the 
former case 7 > 0, the width of 1/3 plateaux has a different 
relation with couplings 7 and 7'. By identifying 6e = ~7' 
and 8e, - —\h, h c = 47' is obtained, that is independent of 
7. Here, the spin configuration on m = 1/3 plateau is illus- 
trated in Fig. [T](c). The order parameter characterizing this 
phase is the spontaneous magnetization m\f,=o, which implies 
the breaking of Z% symmetry. By increasing temperature, the 
Z2 symmetry will eventually recover above the critical temper- 
ature T c , and the spontaneous magnetization will vanish. As 
shown in the inset of Fig. 0(b), the magnetization at T > T c 



starts from m = 0, and the 1/3 plateau is smeared and finally 
destroyed by strong thermal fluctuations. 

In order to look at the effects of external magnetic fields, 
the phase diagrams at zero temperature are plotted in Figs. 
0c) and (d). For 7 > 0, as shown in Fig. 0(c), there are 
four different phases including ferromagnetic, ferrimagnetic, 
plateau ferrimagnetic, and disorder phase that only exists in 
the h — line. It is worthwhile emphasizing that although 
the magnetization in the 1 /3 plateau ferrimagnetic phase has 
the same value as that in the ferrimagnetic phase at T = 0, 
they are quite different in nature. The former is induced by a 
magnetic field and highly degenerate with the degeneracy K = 
2 Nc , while the latter is an spontaneously ordered phase with 
the Z2 symmetry breaking. As indicated in Fig. 0(d), for J < 
0, only a ferrimagnetic phase and a ferromagnetic phase exist. 
Here we would like to stress that the phase transitions between 
these different phases only occur at T — 0, and the temperature 
would then blur the transitions. In fact, in the presence of 
an external magnetic field, there are no thermodynamic phase 
transitions at finite temperatures, which will be discussed in 
Sec. VI. 



B. Susceptibility 



In order to understand the magnetic response of the present 
system to an external magnetic field, the zero-field suscepti- 
bility x 1S obtained by x — \m(K) - m(h = 0)] /h for a small 
magnetic field. In Fig. |6](a), where 7=1 and J' e (0, 1), 
the ground state is disordered, and the spins on one (a) sub- 
lattice are free to flip up or down without an energy cost due 
to the frustration effect, x diverges, obeying Curie law, i.e., 
X °c l/T as T approaches zero. This result is validated in 
the inset of Fig. [6] (a), where the xT curves converge to a 
constant 1/3 at low temperatures, which is independent of J'. 
The specific value 1 /3 can be attributed to the free spins on 
one of three sublattices. On the other hand, in the high tem- 
perature limit, x decays with the Curie- Weiss law, which can 

be fitted by xT = jTS U P t0 T I J = 100 ( Note that Fi 8- 
shows only to T/J 9). It is straightforward to use the mean- 
field approximation to obtain the Curie- Weiss temperature 6 
as (87' + 47)/3, and 6 = 1.867 and 2.933 for 7' = 0.2 and 0.6, 
respectively. The fittings in the inset of Fig. |6](a) agree with 
the mean-field predictions, that further validates our TRG re- 
sults. In Fig. 0(a), it is interesting to notice that there exists a 
turning point at an intermediate temperature in the crossover 
region, which separates the low T Curie behavior and high 
T Curie-Weiss behavior, as shown in the inset of Fig. |6ja). 
Quite differently, as seen in Fig. |6](b), when the ground state 
is ferrimagnetically or ferromagnetic ally ordered, the suscep- 
tibility has a divergent peak at T c where the magnetic ordering 
is destroyed by thermal fluctuations. This again certifies the 
existence of magnetic order-disorder phase transitions. 
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FIG. 7: (Color online) A comparison of TRG results to the experi- 
ment, where the experimental data are taken from Ref. 9i The TRG 
result is calculated at a small field h/J = 0.05. The high tempera- 
ture fittings with the Curie-Weiss law to both experimental and TRG 
results are also shown. The inset predicts a divergent peak in the 
specific heat around T = 20K. 



C. Comparison to Experiments 

The complex Co(N 3 ) 2 (bpg> DMF 4/3 reported in Ref. |g| is 
a molecular magnetic material, in which the Co +2 ions form 
a distorted Kagome layer. Experimentally, the susceptibility 
does not go to zero as T approaches zero (see Fig. |7J, which is 
unusual for an isotropic Heisenberg antiferromagnetic system. 
In fact, the Co +2 ions are believed to have effective spin- 1/2 
when T < 20K, with anisotropic Lande g factors ( g\\ + 0, 
g ± ~ 0), which implies that in this compound the Ising-type 
couplings may be dominant between Co +2 ionsi^ 2 ^ Here, we 
try to use our TRG results to fit the experimental data of sus- 
ceptibility (especially for the low T region) for this complex. 
To be consistent with the experimental convention, the defi- 
nition of susceptibility x = m(h)/h is adopted. As shown in 
Fig- El X decreases steeply around the transition temperature 
and, one may see that the fittings agree rather well with the ex- 
perimental data at low temperatures. The exchange coupling 
constants for this compound are estimated through the fittings 
as J — 22K and J' = 33K. According to our study on the 
Ising DK lattice with the parameters J > and J'/J = 1.5, 
the system has a ferrimagnetic phase at low temperatures. It is 
thus not difficult to understand why the low temperature sus- 
ceptibility goes to a finite value instead of zero for this com- 
pound. At high temperatures, by fitting the TRG results with 
the Curie-Weiss law^ = C'/(T + 6), we find that the Curie- 
Weiss temperature 6 a 161. 3K, which agrees well with that 
of experimental estimation (6 « 165.8/T, see online support- 
ing material of Ref. |9|). Besides, we find that the ratio of the 
experimental susceptibility to the result from the Ising model 
equals a constant R » 5.4 in the high temperature limit. This 
constant ratio may be ascribed to the fact that in the material 
at T >20K the effective spin of Co +2 ions may no longer be 
1/2 and also, the other interactions such as XY couplings may 
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FIG. 8: (Color online) The specific heat C and the magnetization m 
as functions of temperature in different magnetic fields, (a), (b) and 
(e) illustrate the ferrimagnetic case J > 0,J'/J = 1.5, (c) and (f) 
are for the paramagnetic case J > 0, J'/J = 0.5, and (d) depicts the 
ferromagnetic case J > 0, J' / J = -1.5. 



intervene, giving rise to that the Ising model is insufficient to 
describe the behaviors of this complex. Surely, more experi- 
mental results towards this direction are needed. In addition, 
we have calculated the specific heat based on the Ising model 
with the couplings given above, and found that a divergent 
peak exists around T = 20K, as depicted in the inset of Fig. [7] 
suggesting that this compound may undergo a phase transition 
at low temperature. Experimental studies on the specific heat 
and other quantities for this compound will be carried out in 
near future. 



VI. SPECIFIC HEAT AND MAGNETIZATION IN A 
MAGNETIC FIELD 

Next, we will study the effect of an external magnetic field 
on the specific heat. Three typical cases will be studied: fer- 
rimagnetic (J > 0,J'/J = 1.5), paramagnetic (J > 0,J'/J = 
-0.5), and ferromagnetic (J > 0, J'/J = -1.5) cases. 

In Fig. [8] the specific heat in the presence of an exter- 
nal field for the ferrimagnetic case (J > 0, J'/J > 1) is 
shown. Fig. [SJa) shows at small fields with h < 27, the 
peak of specific heat moves towards high temperatures, and 
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its height firstly decreases, and then increases with enhanc- 
ing the field until it approaches the spin flop critical field h c , 
which polarizes all spins. In Fig. 13b), when the field keeps 
increasing, the peak of specific heat becomes dulled, and then 
splits into double peaks (except for the point h = h c ), which 
can be viewed as a field-induced splitting. Similar phenom- 
ena have also been observed in other Ising and Heisenberg 
spin systemsi2i2£ The double peak scenario will eventually 
be spoiled by further increasing h. When h » h c , the spe- 
cific heat will again be single-peaked. It is notable that the 
divergent peaks at zero field disappear immediately when the 
field is switched on, which means that the phase transitions 
are absent and the system remains in the ferrimagnetic phase 
at all temperatures. When h — 0, the ferrimagnetic ordered 
phase spontaneously breaks the Z2 symmetry contained in the 
Hamiltonian (see Eq. []]), and possesses a nonzero order pa- 
rameter. When T > T c , the thermal fluctuations will destroy 
the magnetic order, while Z2 symmetry will be recovered, and 
m vanishes immediately [the solid line in Fig. [8](e)]. How- 
ever, the external field explicitly breaks the Z2 symmetry in 
the Hamiltonian, and m is nonzero even at high temperature 
T > T c [the symbol lines in Fig. [8](e)]. Therefore, no phase 
transition occurs in the presence of a magnetic field. In Fig. 
[8te), according to the magnetization m at zero temperature, 
the curves can be classified into three classes. When h < h c , 
the curves start from m = 1/3; while h > h c , the spins are 
polarized and m = 1 at T = 0; when h = h c , the case is of 
a little subtlety, where m equals to the statistical mean value 
2/3 at zero temperature. 

In Figs. |8](c) and (f), the paramagnetic case J > 0, J' /J - 
0.5 is studied. The field will firstly promote the peak height of 
the specific heat, and moves the peak to the high temperature 
side. When h is close to the critical field, the specific heat will 
again be dulled, where the height is decreasing, and the peak 
splits into two sub peaks except at the point h = h c . When 
the field h » h c a single peak of the specific heat recurs. The 
m—T curves in Fig. [8](f) are quite similar to those in Fig. [8](e) 
and can be classified analogously. At last, the ferromagnetic 
case J > 0,J'/J < -1 is shown in Fig. 0(d). The divergent 
peak for the ferromagnetic-paramagnetic phase transition dis- 
appears owing to the same reason in the ferrimagnetic case 
as mentioned above. The specific heat reveals a round peak, 
which moves towards the high temperature side. By contin- 
uously enhancing the field, the height of the peak decreases 
down firstly and then goes slowly up. 

Besides, we have also studied other situations with differ- 
ent couplings, and found that they can be ascribed into the 
above three classes. For instance, other ferrimagnetic cases 
with J < 0,J' > and ferromagnetic cases with all ferro- 
magnetic couplings (J, J' < 0) behave similarly with those 
presented in Fig. [8] 



VII. SUMMARY AND DISCUSSION 

In this article, we have systematically studied the thermo- 
dynamics and magnetic properties of Ising model on a DK 
lattice by exact solutions and the TRG numerical method. 
It is shown that the phase diagrams are composed of three 
phases including ferromagnetic, ferrimagnetic, and paramag- 
netic phases. Phase transitions between them are identified by 
studying the specific heat and magnetization. The critical ex- 
ponent 77 of m near T c is determined as 1/8, which appears to 
fall into the universality of the 2D Ising models. The TRG re- 
sults of zero-field specific heat agree very well with the exact 
solutions, showing that TRG is an efficient and accurate tool 
in dealing with 2D Ising models. The TRG method is also 
utilized to study the properties in the presence of a magnetic 
field. In the magnetization curves, 1/3 plateaux at low T are 
identified and, the relations of the plateau width with coupling 
constants J, J' are obtained. In addition, the zero tempera- 
ture /' - h phase diagrams are presented to clarify the various 
ground state phases in external magnetic field. The zero-field 
susceptibility ^ of the paramagnetic case (J > Q,\J'/J\ < 1) 
is found to obey Curie law at low T and Curie-Weiss law at 
high T. While in the ferrimagnetic or ferromagnetic case, 
the divergent peak of x is found at the critical temperature. 
Moreover, the specific heat under different magnetic fields is 
also investigated. It is uncovered that the phase transitions 
are absent immediately when a magnetic field is switched on, 
and the field-induced peak splitting of the specific heat is rec- 
ognized when h is close to the critical field. We have also 
fitted the experimental data of susceptibility of the complex 
Co(N 3 ) 2 (bpg)- DMF4/3 with the TRG results, and obtained 
the couplings J = 22K and J' - 33K. Based on TRG calcu- 
lations, a ferrimagnetic-paramagnetic phase transition is ex- 
pected to occur at about T = 20K in this complex, which will 
be studied in future. The present study offers a systematic un- 
derstanding for physical properties of the 2D Ising model on 
the DK lattice, and will be useful for analyzing future exper- 
imental observations in related magnetic materials with DK 
lattices. 
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